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ABSTRACT 



We advocate a new approach to the study of unitary matrix models in external 
fields which emphasizes their relationship to Generalized Kontsevich Models (GKM) 
with non-polynomial potentials. For example, we show that the partition function 
of the Brezin-Gross-Witten Model (BGWM), which is defined as an integral over 
unitary N x N matrices, J[<it7]e Tr ( jtu+Jut - ) , can also be considered as a GKM with 
potential V(X) = j?- Moreover, it can be interpreted as the generating functional 
for correlators in the Penner model. The strong and weak coupling phases of the 
BGWM are identified with the "character" (weak coupling) and " Kontsevich" (strong 
coupling) phases of the GKM, respectively. This sort of GKM deserves classification 
as p = —2 one (i.e. c = 28 or c = —2) when in the Kontsevich phase. This approach 
allows us to further identify the Harish-Chandra-Itzykson-Zuber (IZ) integral with 
a peculiar GKM, which arises in the study of c = 1 theory and, further, with a 
conventional 2-matrix model which is rewritten in Miwa coordinates. Inspired by the 
considered unitary matrix models, some further extensions of the GKM treatment 
which are inspired by the unitary matrix models which we have considered are also 
developed. In particular, as a by-product, a new simple method of fixing the Ward 
identities for matrix models in an external field is presented. 
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1 Introduction 



The full partition function, i.e. the generating functional of all of the relevant correlators of 
an eigenvalue matrix model, can be reduced to an integral over matrix eigenvalues only. In 
this process, the integration over "angular variables" is simultaneously factored out of all of 
the correlation functions. Such models appear to exhibit a novel integrable structure, which 
is intimately related to certain topological structures and which has been the subject of much 
investigation over the last few years. There are essentially two interesting types of eigenvalue 
models which are in fact closely related: one is the family of Hermitean multi-matrix models 
and the other is the generalized Kontsevich Model (GKM). 

However, not all of the interesting matrix models are of the eigenvalue type. For example, in 
models which describe lattice gauge theory, the angular degrees of freedom of unitary matrices 
represent the gauge (Yang-Mills) bosons. Such degrees of freedom are absent in an eigenvalue 
model. Moreover, in the context of string theory, it is thought that non-eigenvalue models should 
be studied in order to overcome the d = 2 ("c = 1") barrier which separates purely topological 
theories from those with non-trivial perturbative spectra. A very interesting and important 
direction for the investigation of these models is to examine in what sense the integrable structure 
which is important in understanding eigenvalue models is still present, if at all, in non-eigenvalue 
models. The full answer to this question, which in the future should probably involve a "non- 
Cartanian" generalization of today's concept of integrable hierarchies of the (multicomponent) 
KP and Toda type, is as yet unknown. (See ref.[l] for more detailed presentation of the problem.) 

A broad and important class of non-eigenvalue unitary matrix models occurs in the situation 
where one is interested in correlators of arbitrary matrix elements of the unitary matrices. 
In that case, an action which depends only on the eigenvalues of the unitary matrices is not 
sufficient to form a generating functional for all of the correlation functions. This sort of model 
is essentially of the non-eigenvalue type. Conventional analyses of unitary matrix integrals 
usually lead to representations in terms of group theoretical quantities (such as characters and 
their multi-variable analogues) which, although intimately related to integrability theory, are 
not transparent enough to reveal the relevant structures. Results of this kind are particularly 
unsuitable for taking large-iV limits, which are often important for relating matrix models to 
physical reality. 

Another important feature of unitary matrix integrals, familiar from studies of lattice Yang- 
Mills theories, is the occurrence in the infinite iV limit of a third order phase transition separating 
weak and strong coupling phases. As first emphasized in ref.[2] two analogous "phases" are in 
fact also present even at finite N in the GKM. The phase structure of unitary matrix models 
thus complements, rather than contradicts their relationship with the GKM and integrability. 

In this paper we describe what we believe is a natural step toward understanding integrable 
properties of non-eigenvalue models. We investigate the simplest example: the Brezin-Gross- 
Witten Model (BGWM) [3, 4] with partition function 

^bgwm(J, J f ) - ±- [ [dU]e^ u+Ju] \ (1.1) 

where the integration is over N x N unitary matrices with Haar measure [dU] and Vjv = 
J NxN [dU] is the volume of unitary group. The "external field" J is an N x N complex matrix. 
This model is a priori of the non-eigenvalue type, since there is no way to reduce the integral 
in (1.1) to that over eigenvalues of U, and also because derivatives with respect to J and jt 
generate correlators of any matrix elements of U and W. Still this particular model is simple 
enough to fit into the traditional integrability framework. Because of the invariance of Haar 
measure Z(J, jt) depends only on N rather than N 2 parameters: the eigenvalues of the matrix 
M = JJ^. Thus, it is not a complete surprise that it can also be represented as the partition 



function of an eigenvalue model, namely of a certain example of a GKM: 



Zn{0) (1.2) 



Z N ( M ) = -1- / dXe Tr(MX-7VlogX + i) 
VjV J NxN 



The integral now is over NxN Hermitean matrices X. 

Being a particular example of a GKM, this model satisfies different kinds of Virasoro con- 
straints in the weak and strong- field limits which, in the terminology of ref.[2]), represent the 
"character" and "Kontsevich" phases, respectively. It is integrable in both limits and, when 
expressed in terms of the appropriate variables (t£ = j-TrM k or t k ^ = |TrM _fe / 2 ) Zbgwm it 
is the r-function of the Toda chain. In the Kontsevich (strong-field) phase, which corresponds 
to the most interesting weak-coupling phase of the original BGWM, it is in fact a reduced - 
MKdV r-function, which was studied earlier in [5], while the relevant Virasoro constraints were 
first derived in [6]. 

This model has various representations which display its connection to other theories. It can 
also be rewritten as 

Z N (M) = / dXe ^x-Ni og x+j±) (L3) 

or, after the change of variables X — > Y = 1/X, as 

Zn{M ) = ^-[ dY e T < M Y- N ^ Y+Y \ (1.4) 
Vn Jnxn 

Representation (1.3) is in fact nothing but the original integral (1.1) with U substituted by 
X and the Haar measure [dU] - by the left- and right- invariant measure 1 



dX 

(Ditxy 



(dX) ^tk^W- (1-6) 



However, this "change of notation" (or of integration contour) reveals an association with the 
GKM and allows one to put the consideration of the BGWM, - and, consequently, of generic 
unitary matrix integrals - into a context which is well adjusted for the study of integrable 
structures. 

The representation (1.4) is particularly important from this point of view. In this represen- 
tation it is possible to "untie" the size of the matrix Y from the size of the matrices J and J', 
substituting (1.4) by 

Z +{L) = ±-\ d Ye^ L ^ Nl ^ Y+Y \ (1.7) 

If it is considered as a function of the t£ = ^tvL k , this quantity is essentially independent of n 
(though it depends non-trivially on N), and thus can be used to reproduce Z^{M). Occurrence 

lr rhis measure has the properties: 

(d(GX)} = (d(XG)) = (dX) for any G; (d-^) = (dX), (1.5) 

X 

which are characteristic of a Haar measure, thus the only real difference (if at all) between X and U integrals can 
be the reality condition (the choice of integration contour). Furthermore, to avoid any confusion, we remind the 
reader that in the theory of the GKM the word "Hermitean" does not mean much more than the "flatness" of 
the measure dX = Y[ { ^ dXij . In many cases it is assumed that the integrals go along some appropriate contours 
in the complex space (of matrix eigenvalues) rather than over real hyp ersurf aces. The matrices can be complex 
in this sense, they should be diagonalizable by unitary conjugation but the eigenvalues are allowed to be complex 
in the way appropriate to making the integral in the GKM well defined. 



of the extra parameter n is important for establishing relation both to the GKM and to Penner 
theory [7, 8]. 

Though it is defined as a single matrix (one-link) integral, the BGWM is known to have 
much to do with 2d Yang-Mills theory [3]. It can be also relevant to d > 2 theories, if considered 
in the mean-field approximation. An example of more sophisticated unitary matrix model, is 
provided by the Harish-Chandra-Itzykson-Zuber (IZ) integral. Though still much simpler than 
generic unitary matrix actions, which are usually defined on the plaquettes, the IZ theory can 
be used for construction of some oversimplified gauge theories beyond 2 dimensions. Making 
use of our results for the BGWM, this integral can be rewritten as a GKM with potential of 
generic shape and, further, as an Hermitean 2-matrix model: 

V N JnxN \ Ui<j Wi ~ <Pj ) Ui<j {<Pi -<Pj) J 

Vn JNxN Vn JNxN 

= c N f [ d HdHe™ R e T ^>-^ HTk+f * nTk) 

J J NxN 

(1.8) 

The two sets of times (with "plus" or "minus" superscript), which appear at the r.h.s. 

T± = iTr$ ±fc , f± = ^-Trl> ±fe , k > 0, (1.9) 

are relevant for description of the two different phases and are in a good sense conjugate to 
each other. 2 Relation (1.8) can be used in both directions: to get a reformulation of the IZ 
integral in terms of a 2-matrix model or to describe the 2-matrix model in terms of the Miwa 
parametrization. This relation reveals the interpretation of the IZ integral as a Toda-lattice 
r-function. 

Since our main goal in this paper is to express the BGWM in terms of the GKM, we begin 
in Section 2 with an exposition of the general theory of the GKM. We emphasize the issues 
of different limits (phases) of the model as well as the peculiarities of the polynomial and 
"antipolynomial" potentials. The latter have not been dealt with in previous investigations of 
the GKM. 

In Section 3 we return to the BGWM and give the proof of the above-formulated statements 
describing its connection to the GKM. Being a very simple example of an antipolynomial GKM 
(the ll p = —2 model") it can be analyzed in somewhat more details than the generic GKM 
which is discussed in Section 2. The issue of Ward identities is particularly straightforward. In 
the p = — 2 case they reduce to Virasoro subalgebras in both the character and the Kontsevich 
phases. 

Finally in Section 4 we present some preliminary discussion of the IZ theory from the point 
of view of the GKM. Further development in this direction should involve consideration of cor- 
relators with the IZ measure and leads beyond the standard mathematical notion of "cartanian 
integrable hierarchies" . 

Our conclusion emphasizes, among other things, that formalism of this paper allows us to 
identify the "phases" of the GKM with those of physical gauge theories. 

2 Definition of the "zero-times" deserves a comment: actually the item with k = in the sums in the 
exponent in (1.8) should be ascribed different meanings for different choices of superscripts: 

T+H~° -> -Tr(J (g> \ogH) = -A^Trlogif; 
while TqH +0 -> -Tr(log $ I) = -NTr log 



2 Two "phases" and two versions of the GKM 



Let us recall that the main object in the theory of the GKM is the "Kontsevich integral" 

Tv{Ni L) = * f dXe tr(LX-mo g X + V(X)) (21) 
V n Jnxn 

over n x n Hermitean matrices X with an Hermitean external field L and some "potential" 
function V{x). This is actually an eigenvalue model, since both matrices L and X can be diag- 
onalized and the angular variables can be integrated away with the help of the IZ formula (first 
line of (1-8)). We refer the reader to [1] for a review of the basic theory of Kontsevich integrals 
and a list of the relevant references. For the purposes of this paper some novel considerations are 
necessary. In particular, it is important to distinguish between four different situations, where 
the issues of integrability and Ward identities deserve separate analyses. 

First of all, it is important to know where the singularities x s of the potential lie. We shall 
distinguish two cases: x s = oo (the "polynomial" case) and x s = (the "antipolynomial" case). 
Actually in this paper we mostly deal with monomial potentials, V{x) = — where p can 
be either a positive or a negative integer. Below we use P to denote the absolute value of p: 
P=\p\. 

Another distinction is between two different limits ("phases") of the GKM. These can be 
viewed simply as two different asymptotics of the Kontsevich integral (2.1): for small and large 
external fields L. The term "phases" is not really accurate in the theory of the GKM where the 
infinite-dimensional phase space is analyzed. This terminology derives from the fact that the 
weak- and strong-field limits of (certain versions of) the GKM are identified with the strong- 
and weak-coupling "phases" of lattice Yang-Mills models. 

The currently-available results about the four subtopics of the GKM theory are summarized 
in the following table: 



Table 1 





polynomial 


antipolynomial 


character 
phase 


? 

yy ("^-constraints 

? 


? 

y\? (+, \p 1 ) -constraints 

? 


Kontsevich 
phase 


p-reduced KP r-function 

(r) 

Wp k -const., k > 1 — r, r < p 
string eq. C- p t = 


\p -reduced KP r-function 

(r) 

Wp k -const., k > 2 — r, r < \p\ 
np\) 

string eq. %( p | 2) r = 



Only one of the four cases has been thoroughly investigates, that of the Kontsevich phase 
of the polynomial model. The structure of the character phase has so far been addressed only 
in ref.[2], but knowledge of its detailed structure is thus far incomplete. For comparison with 
the BGWM we require, for the most part, the antipolynomial model which to the best of our 
knowledge has never been analyzed in detail. We shall find that, especially in the case of the 
antipolynomial potential V(x) = l/x which is most relevant to the BGWM, some results are 
accessible. The purpose of this section is to describe the content of table 1 in a little more detail. 
Many things remain obscure and deserve further investigation. 

For the BGWM only the piece of a generic GKM theory which deals with monomial potentials 
is actually relevant, therefore we often restrict consideration in this section to monomials, thus 
leaving a very interesting piece of the GKM theory beyond the scope of this paper. 



2.1 Time- variables 



Kontsevich integral J-y(L), introduced in (2.1), is a symmetric (Weyl-group invariant) function 
of eigenvalues l a and has different asymptotics for small and large Vs. In this paper we discuss 
only the two simplest situations: when all the Vs are either large or small at once. In the first 
limit ("character phase") J-'v(L) can be considered as a function of Weyl-invariant quantities 
t+ = ±trL fc , and 

Z£ KM (t + \n,Af,V) = T V {M,L). (2.2) 

In the second limit ( "Kontsevich phase" ) the integral is first rewritten in terms of eigenvalues 
of X and L and then is evaluated by the saddle-point method. The saddle-point xq(1) is defined 
from the equation 

l + V(x (l)) = 0. (2.3) 

For large I the eigenvalues xq(1) are close to singularities of V'(x) and much depends on the 
degree of the singularity. If V'(x) ~ (x — x s )~ p , then xo — x s ~ l~ l l p and perturbation theory 
around the saddle point is expandable in terms of 3 

In order to define ZQ KM (t~) it remains to factor out the quasi-classical contribution to ^(L): 
Z GKu(t- {p) W,V)\ t{ _ p)= ^ tTL _ k _, = (Cv(AA,L))- 1 ^ v (AA,L), (2.5) 

kp+i k + i/p 



where 4 



e tr(LXo+V(X )) 



2.2 Integrable structure 

Integrable structure (of cartanian type) of the GKM is essentially due to existence of determinant 
representation 

*w. L » = A© n / 1^* 1 ^") = ""-y^ (2 - 7) 

where A(x) = Ila>b( a; a - x b ) and 

* (Z) = J x a - x e xl+v{ ^dx. (2.8) 

This expression should be compared with the standard representation of (1-component) KP 
r- function in Miwa coordinates, 



r , -i det ab [\fy a -f3(\b)} 
T/3\tki = — 



A(A) 



(2-9) 

tk ~ k Ea=l ^ a 



We reserve a simplified notation t k for t k , i.e. for p = 1. Everywhere below fc > 0, but the zero-times to are 
defined in a somewhat tricky way: ~ ]jtr/; ~ h^ x L~ k as fc — > +0. This means that £q is usually substituted 
by t ( y = — trlogL, and do not contribute to the sums ~^2 k>0 kt k , while t£ on the contrary, does contribute, and 
kt^ \ k =o = trl = n. 

4 Note that logarithmic piece is not included into potential V(x) and is treated separately in (2.3) and (2.6). 
This prescription is important for the nice integrability properties of Z^ KM . Note that it is self-consistent provided 

Af/xo~l for large I (what will usually be the case in our considerations below). 



As usual, this formula describes the restriction of r-function on the n-dimensional hypersurface 
in the infinite-dimensional space of time- variables, and this is the only source of n-dependence 
at the l.h.s. as the shape of the function r{tk} does not depend on any parameter like n. r- 
function depends on choice of the (n-independent) functions (called basis vectors) ip, which are 
only restricted to have definite asymptotics at large A: 



^ a (A)~A a - 1 (l + 0(l/A)), (2.10) 

where the r.h.s. should be expandable in negative integer powers of A. Classes of equivalence of 
such ip (modulo analytic at A = oo transformations of A) label the points of infinite Grassmannian 
(the "universal module space" ) . Parameter f3 is referred to as " zero-time" . It can be considered 
as a remnant of embedding of KP hierarchy into a more general Toda-lattice hierarchy. The 
hierarchy is P-reduced, whenever 

a+P-l 

A P Va(A) = Va+P(A) + Yl Q<*MX) i 2 - 11 ) 

6=1 

with any A-independent Q a b. Characteristic feature of the P-reduced r-function is that 



~ T ^ ~ = cons ~ti f° r an y integer k,l, (2-12) 

i.e. the r-function is essentially independent of all the time-variables tp^- 

In order to obtain the r-function interpretation of the GKM it is only necessary to adjust 
the ^-functions (2.8) so that they satisfy conditions (2.10). For this purpose one can connect 
properly I and A, change normalization of functions \l/ a and also change the labeling (a indices). 
Besides, we should do linear combinations of 4> a 's not changing the determinant (3.30). There- 
fore, we can add to the second vector the first one with an arbitrary coefficient, to the third 
one both the first and the second ones with arbitrary coefficients and so on (so it corresponds 
to triangle substitution which does not change the determinant). Then we adjust these linear 
combinations such that new vectors ^ a (l) will have the asymptotics l a ~ l . Let us also normalize 
them to have the unit coefficient in the first term. 

To conclude these introductory comments, let us briefly explain why the determinant repre- 
sentation (3.30) with V>a(A) independent of n and the asymptotics (2.10) guarantees n-independence 
of (3.30) as a function of times. In fact, it means that changing n + 1 — > n does not change the 
functional form of (3.30), but only the point in the time space, i.e. times are now parameterized 
by n instead of n + 1 Miwa coordinates. Put differently, this means that one can tend A ra+ i to 
oo, and the determinant in (3.30) will be of the same form, but ofnxn matrix. This is actually 
the case provided by the condition (2.10) and unit coefficient in the leading term. This unit 
coefficient is also consistent with the normalization of (2.1) which implied to be unit when all 
A tends to oo, i.e. all the times are equal to zero. We postpone more detailed discussion to the 
conclusion of the paper, when we will have some manifest examples of how it works. 



2.2.1 Character phase 

In the character phase ^ a (l) are expanded in positive integer powers of I, therefore A should be 
identified with l/l. However, * (Z) as defined in (2.8) is * a (Z) ~ 1 + 0(1 = 1/A) rather than 
the what is required in (2.10). However, we are still free to perform a linear transformation of 
the set {^ a } which does not change the determinant (3.30). Namely, we can add to the second 
vector the first one with an arbitrary coefficient, to the third one both the first and the second 
ones with arbitrary coefficients and so on (thus the linear transformation is triangular and leave 
determinant intact). These linear combinations can be adjusted so that the new vectors 4> a (l) 



have asymptotics l a . Then, we can also normalize them to have unit coefficients in front of 
the leading terms: this gives rise to an overall A-independent constant. 

For these new functions we have: 

i! a . N {l) = l a - M -\l + 0(l)), (2.13) 

but this is still "inverse" with respect to (2.10): A 1_a appears instead of A a_1 . Now comes 

the crucial step, which one should always make in the character phase [2]: let us relabel the 
^-functions according to the rule 

§ a = * n _a+i. (2.14) 

i.e. just replace a — ► n — a + 1. This transformation does not change the determinant: 

det * a (Z 6 ) = det * a (Z 6 ) = det l a (Z b ) (2.15) 

ab ao ah 

Then 

l a = l n ~ a (i + 0{l)) = A 1_n [A a_1 (l + 0{\))} = X^MX) (2.16) 

and taking into account the change of Van-der- Monde determinant in (2.9) when I — > A -1 , one 
finally obtains 

j> a = \ a - 1 (l + 0(\)). (2.17) 

These functions already possess correct asymptotics behavior (2.10), but instead they can depend 
on n. Indeed, starting from n- independent functions we performed the transformation (2.14) 
which introduced manifest dependence of n. In general situation this is a disastrous obstacle 
for the character phase partition function to be r-function. Nevertheless, sometimes by special 
adjustment of the coefficients in the potential in (2.8) it is possible to get rid of this n-dependence. 
An example of such a situation is provided by BGWM and will be considered in section 3.2.3 . 

5 



2.2.2 Kontsevich phase. Polynomial model 

In the Kontsevich phase things will be somewhat different in the polynomial and in the an- 
tipolynomial situations. Let us begin with the polynomial case which was considered in detail 
in [9, 10, 11]. The Kontsevich phase corresponds to a large I expansion. Therefore, using the 
saddle-point method, the integral (2.8) can be expanded in a series of inverse powers of I. Indeed, 
it is straightforward to derive 

5 Within the character phase it is also interesting to consider discrete sums (over eigenvalues of X) rather 
than integral in the definition of TviL). Such discretized GKM is intimately related to representation theory of 
compact groups and in this context it can be also reasonable to consider expansion in powers of f£ = |;tre fei = 

fi-iy. tf • ^ ' n Kontsevich phase the saddle point occurs at the large value of x — xo, then there is no real 
difference between discretized and continuous GKM. While this condition is satisfied for potentials V(x) which 
are polynomials in x (so that x s — oo and x (l) ~ l +1/p ), it is no longer true for polynomials in x 1 (i.e. x s = 
and xo(l) ~ l^ 1 ^), which is actually the case in our discussion of the BGWM: see eq.(3.9). In this case, however, 
the better representation is in terms of Y = l/X, and yo{l) ~ l 1 ^ is large in Kontsevich phase. 

Yet another possibility to represent (2.8) in the form of a r-function was proposed in [2], where all n-dependence 
of basis vectors is ascribed to n additional Miwa variables parameterizing the negative times of Toda lattice 
hierarchy. Then, basis vectors are n-independent but instead only on negative times. However, putting these 
new parameters to be unity, one gets the coefficients of basis vectors manifestly depending on n. In this case 
the n-dependence is completely due to the special choice of the negative times, and the determinant (3.30) still 
describes r-function, but that of Toda lattice hierarchy. 



*a-JV(0 = Jv"(X)e xv '^- v ^\- M X a - 1 (l + O(-)) = 

V A (2.18) 

= V / V^(A)e AV '( A )- v ( A )A-^[A^a-^(A)] ^ ,- 1 (A)[A^ a _ A r(A)], 
where A is solution of the saddle-point equation: 

l + V(X) = 0. (2.19) 

Then, in accordance with our general rules, we should suit normalization of \I/ a to the asymptotics 
conditions (2.10). This normalization factor in matrix form is just Cy of (2.6). Simultaneously, 
we automatically obtain basis vectors which do not depend on n and the coefficient M in front 
of logarithm plays the role of " zero time" . 

Consider the case when the potential is a monomial of degree p + 1: V = — ^+37- Then, by 
examining small variations of x in the integral (2.8), one obtains the recurrence relation 

Z* a (A) - * a+p (A) = (a - 1) a _- (A) (2.20) 

and a similar one for the functions ip a (\). Using (2.19), i.e. I = X p , one immediately gets the 
reduction condition (2.11): 

A p Va(A) = Va+ P (A) + (a - l)Va-i(A). (2.21) 

Thus, we prove that the GKM in the Kontsevich polynomial phase with monomial potential of 
the degree p describes the r-function of p-reduced KP hierarchy. 

We demonstrate now that antipolynomial model is nothing but "analytical continuation" of 
the polynomial model to negative values of p. 



2.2.3 Kontsevich phase. Antipolynomial model 

Thus, let us consider antipolynomial Kontsevich phase, i.e. potential V(x) in (2.8) which is 
polynomial in inverse degrees of x. Again we can apply the saddle-point method to get the 
same saddle-point equation (2.19) which usually allows one to connect I and A, and we can read 
off the correct normalization factor from (2.6). 

Nevertheless, there is one subtlety which requires some accuracy. To demonstrate the point, 

x p+i x ~ p+1 

let us consider again the monomial potential V(x) = -j— j- = — — , where P = \p\ = —p. 

Then, I = Xq = Xq p , where xq is the saddle point, and we choose A = Xg 1 (as the Kontsevich 
phase requires A to increase with I) and we get for the asymptotics expansion of correctly 
normalized functions * a (A) (these are still not correct ijj a (\)): 

§ a _^(A) = A 1 -*e-^ AP "M/ a _ A r(A) = A^- a (l + 0(1/A)). (2.22) 

For given n (the size of determinant in (2.9) this set of functions can be transformed to that 
with increasing power of A, as required in (2.10), by relabelling. Namely, one can substitute 
a — > n — a + 1 and get 

*„- -jv+i(A) = A a - 1+A/ "-(l + 0(1/A)). (2.23) 

These vectors have correct asymptotics behavior, but depend manifestly on n. Therefore, as the 
last step one should choose N = n — (3 to get rid of any n dependence: 



ip a -p(\) = *^- a +i(A) = s(X) J X P- a e xXP+v ^dx. (2.24) 

This actually means that, while the integrand (2.1) contains manifest n dependence, the integral 
as a function of times does not! Simultaneously we get interpretation of the difference n — N as 
the zero time. 

Now let us look at the recurrence relation satisfied by the basis vectors ^> (A). It can be 
derived similarly to (2.20) and looks like 

A p § a (A) - § a _ P (A) = (a - l)* a _i(A). (2.25) 
After relabelling, one obtains the reduction condition (2.11) with P = —\p\: 

A P Va(A) = Va+p(A) - a^o+i(A). (2.26) 

Therefore, we proved that the Kontsevich antipolynomial phase is described by the r-function 
of \p\ -reduced KP hierarchy. This implies that any GKM in Kontsevich phase can be described 
as the r-function of |p|-reduced KP hierarchy, p being the degree of the potential. 

2.3 String equation 

Besides integrabilty, another important property which characterizes the GKM is that its parti- 
tion function satisfies some constraint algebra, which determines it unambigously. In this sense, 
one does not need any specific information about integrability. However, it is usually a very 
difficult problem to find this constraint algebra . Considerably easier is to use integrabilty prop- 
erties (which fix the partition function as a r-function) and one more equation, namely, "string 
equation" 6 , which specifies r-function uniquely (i.e. fixes the point of the Grassmannian unam- 
bigously). The string equation is usually the first constraint of the algebra. We will see in the 
next section that the constraint algebra can be easily obtained in the character phase. However, 
we think it is too early to recognize the first constraint of this algebra as string equation, since to 
acquire its full status it should be sufficient for restoration of the whole constraint algebra when 
supplemented by appropriate integrability property. But these properties are still undiscovered 
in the character phase. This is why there remain empty places in the table 1. 

Unlike this, in the Kontsevich phase we know the integrable properties, but the entire con- 
straint algebra can be found only after some tedious work. Therefore, in this case, the string 
equation, which can be easily found, is a very efficient tool. Below we derive string equations in 
both models in Kontsevich phase. 

2.3.1 Kontsevich polynomial phase 

In this subsection we shall shortly repeat the derivation of string equation in the polynomial case, 
addressing reader to [9] for more details. We shall demonstrate that this derivation crucially 
uses only the following information: 

1) given asymptotics of tp a (l^) ~ z^ -1 ; 

2) the manifest form of normalization factor depending on M in the GKM integral; 

3) the manifest form of linear term in X in exponential in the integrand in (2.1). 

6 In the literature, it is used to call string equation the derivative of the appropriate constraint over the first 
time variable. For the sake of brevity, we usually call string equation the constraint itself. 



The main idea is to consider the following derivative of the GKM partition function r 
det^(A 6 ) _ dets(X b )* a (X b )_ (gee (2 9)> (218)): 



^\T7^irA^T\ (2.27) 



A(A) - A(A) 

Tr i 

. V"{A) dA 

where, as usual, L = V'(A), and rewrite it as 

v 1 dtp dlogr _ v 1 1 dlogr 

^ Q V"{A) dA dt p £j V"{A) AP+i dt p ■ K ■ > 

On the other hand, the derivative (2.27) is equal to two pieces: the first one originates from the 
derivative of factors s(A) and A(A) -1 (here we use the information of point 2) and is equal to: 

ly 1 V"{K)-V"{\ b ) , . 

2 ^ v"(\ a )V"{\ b ) x a - x b [ ■ > 

where potential V(X) generally contains logarithmic term. The remaining second piece can be 
transformed to derivative over t\ essentially using the correct asymptotics of V'a(A) (point 1): 

T ^vk)^ d ^M = k' ogT - (2 - 30) 

Thus, we finally obtain the string equation in the form: 



' ~h dt p + dt 1 + 2^ b V"(X a )V"(X b ) X a -X b 



(2.31) 



with 



Now let us describe some different approach to the string equation proposed in [12]. For the 
sake of simplicity, we consider only monomial potential. 

The idea of approach is to use the integrability of the system, i.e. the fact that its partition 
function is a r-function. This means that, instead of considering operators acting in the spce of 
functions of times, one can immediately operate in the infinite-dimesional Grassmannian. Our 
problem now is to find the analog of the string equation in these terms to fix unambigously the 
point of the Grassmannian. We already specified the r-function of polynomial model to be a 
p-reduced r-function. This means that the element of the Grassmannian we are looking for lies 
in the subspace V satisfying the condition: 

X P V = IV C V. (2.33) 

Now we are searching for one more condition on V. For this purpose, let us note that all 
functions ^f a (X) are connected by the transformation 



or, equivalently, 



V a+1 (l) = -V a (l), 



(2.34) 



Ml) = AMI)- (2-35) 

The set of basis vectors V'a(A) gives an element of the Grassmannian. On the other hand, the 
operator A maps the set of basis vectors onto itself. Thus, we get the following new condition 
on the subspace of the Grassmannian: 



Mi(i) 



»(() r -(0 



AV C V. (2.36) 

This condition is equivalent to the string equation (see also [12]) and, along with (2.33), specifies 
the poitn of the Grassmannian, i.e. the GKM partition function. In particular, one can see that 
this is actually generating the constraint algebra. Indeed, conditions (2.33) and (2.36) imply 
that any commutator of products of I and A also leaves point of V in V. On the other hand, 
the commutator [A, I] = 1, which means that all possible products l a A b , a, b < generate a 
subalgebra of iy(°°)-algebra. However, this algebra does not annihilate the r-function, as this is 
the algebra in the space of the parameter Since the true spectral parameter is A (i.e. all times 
are made of integer degrees of A), one should make a transformation X p — ► A. It turns out to 
be extremely non-trivial problem [13]. We demonstrate in section 2.5 how to find the constraint 
algebra without explicitly making this transformation. 

To conclude, let us remark that the operator X a ^ (in the space of the spectral parameter) 

corresponds trivially to the operator W^X 1 ^. This means that the operator A corresponds to 
the £_ p -constraint, which is, in fact, the string equation operator A (2.31). Besides, it implies 

that the subalgebra of VT^^-algebra, annihilatig the GKM in polynomial case, contains only 

(r) 

modes which are zero by modulo p, and, moreover, of W^ k -generators, only those with k > 1 — r 
are presented. We return to this issue in section 2.5. 



2.3.2 Kontsevich antipolynomial phase 

Now let us repeat the derivation of the string equation for the antipolynomial model in Kontse- 
vich phase. It is difficult to reproduce the derivation of the first part of the previous subsection 
for the general case of antipolynomial model. For illustrative purpose, we demonstrate this 
derivation in the case of p = —2 in section 3.3.4. But now we are going to get the string 
equation in terms of the Grassmannian. 

In the antipolynomial case we have again the reduction condition 



\ P V = IV C V. (2.37) 

Now we need to construct the operator A. To do this, let us note that the derivative Jj 
shifts the index of the basis vector (2.24) in wrong direction: s(A)^s~ 1 (A) ip a — ► ip a -i, i.e. 
s(X)§ l s- 1 (X) V t V. However, we can use the reduction condition (2.37) to shift it back. 

Thus, the operator / s(A)^s _1 (A) ip a — ► ipa+p~i + where dots stand for terms which can 
be removed by the linear low-triangle transformation of basis vectors. This operator is a good 
candidate for A-operator and it satisfies (2.36). Along with the reduction operator, it generates 
the subalgebra of H^ 00 -* -algebra with leading terms /"J^f with a > b. But this is not the 
maximal possible algebra. Indeed, let us note that we can take as many as P — 1 derivatives 
of basis vector and after this use the reduction condition. This procedure induced by operator 
A = \ls{X)^s-\X) 



operator we need. Sure 



transforms the basis vectors ip a — > ipa+i. This is just the "minimal" 
y, it still satisfies the equation (2.36). 



It mean, finally, that operator A acting on the Grassmannian has the form 



A = 



x -p 2 +3P-i ® P 1 , 
dX P-i i " 



(2.38) 



and we obtain the string equation (2.36). In accordance with the general rules above this 

(P) 

operator corresponds to VF_^p_ 2 ^ p -constraint. It proves the corresponding statement in the 
table 1. 

Let us note that one can again construct the subalgebra of W^°°)-algebra from products of 
the operators I = X p and A: l a A b . It again contains only zero by modulo p modes, but different 

(r) 

restrictions will be produced: in the W^ k -algebra only k > 2 — r modes are presented for r < p 
and k > 2 — p for all other r. Let us note that we can consider a W- algebra defined at the 
vicinity of oo, instead of 0. This will result into the change of sign of all modes. We will use 
in future just this convention (and it is used in the table 1) to have common approach to both 
polynomial and antipolynomial models. 



2.4 Ward identities 

All the Ward identities for the GKM follow from the matrix-valued equation of motion, (L — 
+ V'(X)) = 0, which can be rewritten as 

\l-N (-JL) 1 + V ) ^v(AA, L) = 0. (2.39) 



dL tr J \dL tT 

If M 7^ and/or V'(x) contains some negative powers of x, the reasonable identity arises as 
some L-derivative of this relation, so that it becomes differential rather than integro-differential 
equation (eq.(3.8) is a particular example). 

2.4.1 W-operators 

Ward identities for the GKM in the character and Kontsevich phases arise when ZQ KM (t ± ) from 
(2.2) and (2.5) are substituted for T into (2.39). The resulting equations are expressed in terms 
of the differential W-operators. These are defined by any of the following three relations: 



f> \ m+l 

o£-) /(^) = E £± " 1 w£ ,+1) (W) 

dLto / S >1 



(2.40) 



7 

or' 



or 



w£™ +1 \t)e^>^ LTk = jtr (^) m ^ S } e^>-^ LT \ (2.41) 



n&r 1} W = E kttW^Jt) + £ £^£&n(t)- (2-42) 



fc>0 k=l 
The last recurrence relation should be supplemented by "initial condition" 



or even 



Wi ± ' 1 ) = A, S >1 (2.43) 



W ( s ±fi) = 5s,o. (2-44) 



7 It is in eqs.(2.41) and (2.42) that the convention kt^\k=o = til = n, introduced in the footnote 3 is essential. 



These relations define operators W^'™ (t) for s > 1/2 ± 1/2, no reasonable definition of 
harmonics with s < 1/2 ± 1/2 is known. This and the recurrence relation (2.42) are their 
most striking differences from the conventional Zamolodchikov's operators W(t), defined by the 
standard bosonization procedure. 8 



2.4.2 Character phase 

Relation (2.40) can be used directly to derive the Ward identities in the character phase of the 
polynomial GKM with Af = 0. If V(x) = -j^-, p>0, eq.(2.39) turns into: 

E L 3 ' 1 { {MX P P Ut + ) - 5 S , 2 ) ^km(^)} = 0, ( 2 .45) 

S>1 

from which we conclude, that 

Mtf-i(t + )Z£ KM (t + ) = 5 s , 2 Z+ KM (t+), 8>1. (2.46) 

If potential is not a monomial, a sum over p arises at the l.h.s. 

If Af / 0, one should consider an L-derivative of (2.39) in order to get rid of the integral 
operator (d/dL tr )~ l : 

[A, L - m+ ik/{Jh)}^ L) ^ (2A7) 

This can be again rewritten in terms of the W-operators, but there are two essential differences 
from the case of Af = 0. First, operator VV(+'P +1 ) will appear instead of the W( + < p \ Second, 
explicit dependence on n, i.e. on the size of the matrix L will arise, because when acts on 
L, it produces nl. (It is not just a commutator because the contraction of matrix indices remain 
intact.) Keeping these two remarks in mind, we obtain: 

J2 L*- 1 {(Vtftr V) - W - n)S s ,i + nfo^l - 6, tl )) Z+ KM {t+\N)} = o, (2 . 48) 



S>1 



or 



V) + ^-(1 - Z+ KM (t+\M) = (Af- n)5 M Z+ KM (t+|AA), 

s > 1. 

If Af = this relation is of course a corollary of (2.46) and (2.42). 

Now it is clear, what should be done in the anti-polynomial case. For p < and V(x) = 
— frpT = T^T> P = \p\ = ~P> one should take as many as P derivatives of (2.39): 

' 8 ^-WA)" + mVm) WIH . ( , 50) 



\\dL tr J \dL tr J \dL tT J \dL tr / j 

Then the last item at the l.h.s. is just —I, while the first one can be rewritten as 

d 



a + b=P-l 

a,b>l 



L 



tr 



dL 



tr 



F v (Af,L) = 0. (2.51) 



tr 



8 >V-operators of ref.[14] are in fact the what we now call VV' ''s. The VV^'s have not been discussed in the 
literature so far. 



From (2.40) we deduce that 



E 

a + b=P-l 
a,b>l 



d 



dL tT 



Tr 



dL tI 



-'tr 



a+b=P-l 

a>,b,s ,q> 1 



9 



<9L 



tr 



E^& P i 



tr s>l 

and, putting everything together and using again (2.40), we obtain: 



E^ 1 



S>1 



a + fe=P-l 



or 



+ (n - AOwi^ + W<+£> 2 (1 - - S. tl } Z+ KM (t+\N) = 0, 



E ((« - 1)^1^) + (n - ^)wSf_- 2 1) + 



a +6=P-l 



+ w££ 2 (i - 2+KM(t + Us0 = 5 s ,iZ+ KM (t + |AA). 

Of all the W^-operators in this paper we need only 



s-l 



4 d) (n, t) ee = n J- + £ *t fc ^- + E 



d 2 



9t s k>Q dt k+s k=1 dt k dt s _ k 



(2.52) 



(2.53) 



(2.54) 



(2.55) 



We emphasize, that this operator depends explicitly on the parameter n (the size of L matrix), 
which is hidden in the previous formulas in the k = item in the sums over k (according to our 
convention, ktk\ k=0 = n). 



2.4.3 Kontsevich phase 

We have demonstrated how to generate Ward identities in the character phase and how to trans- 
form them to PF-operators. This is, however, somewhat less straightforward in the Kontsevich 
phase, because one should take into account the Cy factor and the difference between the proper 
variables t( -p ) and t~ = t^~ l \ which appear in (2.40). Indeed, what one should do is to take in 
to account correctly the terms which are constributed Ward identities derivatives of the factor 
Cy. This adds to the standard GKM Ward identities (2.39) some new pieces. This pieces are 
of importance as (2.39) failed to be expanded properly into traces of negative powers of L such 
that the result depends only on times t pk+ i (2.4). But it can be done after taking into account 
all the contributions from Cy! 

Simultaneously it generates the system of VF-constraints imposed on the partition function. 
The drawback of such a calculation is that this is very long and tedious and was presented only 
in the cases of p = 1 [15, 10], p = 2 [16, 17, 6], p = 3 [18] and p = -2 [6]. (We will return to 
the last case in section 3.3.) Thus, we need some other way to determine the constraint algebra. 
One more way was proposed in the paper [13], who manifestly demonstrated the tranformation 
I — ► A at the level of VF-algebra (see discussion in section 2.3). This calculation is also rather 
tedious and was not completed in full. Instead, in the next section we propose a very simple, a 
bit heuristic way to determine the constraint algebra which is imposed on the partition function. 
It allows us, in particular, to determine this in the antipolynomial phase. 



2.5 Ward identities as recursive relations. A new way to deal with Ward 
identities in matrix models 

Now we are going to explain new and to our knowledge the most effective method to deal with 
constraint algebras in matrix models though it is not very rigid, but rather heuristic. 

It is based on a specifics of matrix models that the Ward identities are essentially the same 
as equations of motion and thus define the partition function unambiguously, at least in the 
form of a formal series, 

Z = 1 + ^a s t s + ^2 a si,s 2 ts 1 ts 2 + (2.56) 

S S1,S2 

We shall now explain how this actually works in different models and phases and how it can be 
used to fix constraint algeba. 

All kinds of W-operators (including Virasoro, VV^, Zamolodchikov's W etc) possess the 
following property: operator with the subscript r is a linear combination of terms like 

d b 
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ot h . . . dt h (2.57) 
l\ + . . . + lb — k± — . . . k a = r and all Zi, . . . , lb > 1, 

where i^'s are the times t^s, may be shifted by a constant (see below). A property of the systems 
of constraints, arising in the study of matrix models, which is responsible for the uniqueness 
of their solutions, is that every constrain has an item i.e. with linear derivative and t- 
independent coefficient, and every integer r > 1 appears in one and exactly one constraint. 

The origin of such terms is somewhat different in different situations. 



2.5.1 Character phase 

In the character phase ik = tk, s = s and such terms arise from the contributions with kt\ \ k=o = 

n to VV-operators. appears for this reason in the operator VVr ±,p ^ and preserve the gradation. 
In order to have all the integer r > 1 represented exactly once in the systems of constraints (2.48) 
or (2.53) the labels r in these systems are restricted to be r > 1. 

Whenever such system of constraints is given, it unambigously fixes the perturbative expan- 
sion of the partition function. To prove this, one should just act iteratively starting from the 
first gradation level. It is possible as the coefficients {a Sl) ... )Sfc } has the determined gradation 
level s\ + S2 + ... + Sk and so do constraint algebras (2.48) and (2.53). We illustrate the procedure 
for the simplest example of P = 2 antipolynomial phase. Then we get from (2.54) and (2.55) 

y dt s ^ k dt k+s fr[dt k dt s _ k ) 

As the first step, we calculate a\ using /^-constraint at all t& = 0: 

Na ± = l, ai = ^- (2-59) 

Now we have two coefficients at the second level (in the gradation) - a 2 and a^\. They can be 
fixed by two equations obtained from /^-constraint and by differentiating /^-constraint in the 
first time, both taken at all tk = 0: 



(2.58) 



Na 2 + 2oi,i = 0, 



27Voi,i + a 2 = ai, (2.60) 
1 1 
Le ' 01,1 " 2(JV2-1)' a2 ""iV(Af2-l)- 

This procedure can be evidently continued to build all the expansion (2.56). In particular, at 
the fe-th level there are P(k) unknown coefficients a^, 0^-1,1, Ofc-2,i,i) P(k) being the number 
of partitions of k into integers, which can be determined from P(k) equations: constraints 
]$-£k~i, ■ ■ ■ , all taken at tk = 0. It is clear that the number of unknown coeffiecients 

coincides with the number of equations. 

Besides, it is clear that the same procedure is applicable for any constraint systems of (2.48) 
or (2.53) type (note that no algebras of several different spins simultaneously presented here - 
in constrast to Kontsevich phase below). Indeed, we could guess these constraints even without 
concrete calculations. This is not of great importance in the present trivial case, but gives new 
and very efficient method in the compicated case of Kontsevich phase. 



2.5.2 Kontsevich phase 

In Kontsevich phase there is a "shift of times" ktk = ktk — pSk, p +i, or, better, 

= V ( S - (2-61) 

both for positive and negative p. 9 Such shift breaks the gradation rule (2.57). Nevertheless, 
the reasoning above is still applicable, as this shift goes in "correct" direction, which means that 
it preserves hierarchical structure of equations for the coefficients {o Sl ,..., Sfc } in (2.56) and they 
still can be found iteratively. 

Now let us consider the concrete structure of Kontsevich phase with the deformed grada- 
tion rules given by the shift (2.61). We should take into account that there should appear 

(r) 

only Wp k -operators, as we need the operator algebra which respects p-reduction, i.e. does 
not depend on t, p k and have natural p-gradation. Then, the term arises from the term 

A A f) (r I 1) 

k\tki ■■■Ktu-m - in the W:_,' , ir operator. Actually s — r(p + 1) = pk, thus 

s = r mod p. This is the main complication as compared to the character phase which involves 
separate consideration of the constraints of different spins. Indeed, we need the constraint sys- 
tem with each -J^- term appearing in exactly one equation (for all positive, or for all negative s 
depending on the direction of gradation). 

As the first example, let us consider polynomial case. Then, term can be obtained from 
£_p-constraint, d f +l term - from /^-constraint, ... , and, generally, Q t ® +1 term - from £(k-i)p~ 
constraint. Similarly we can get all other constraints and obtain finally that the constraint 

(r) 

algebra fixing our partition function unambigously is the system of W\ -operators with r = 

l,.-.,p and k > 1 — r. Let us note that the terms with -^-^ arise from the W^ k constraints, 
expressing the independence of the GKM partition function in Kontsevich phase of the t p k- 
variables (this is a little more than just p-reduction - see (2.12)). 

Now let us consider the antipolynomial case. This time dt 9 +1 term again arises from .Co- 
constraint, but we do not need term, i.e. now our /^.-constraints w iH be limited to k > 0. 

9 We discuss here the implications of this shift. As to its origin and relation to topological models and quasi- 
classical hierarchies, we refer to the detailed discussion in ref.[19]. 



Analogously, the term g t +2 appears in VF_ p -constraint and etc. up to W^Jf -constraint which 

is limited this time by k > 2— \p\. It means that this time we have the system of W^-operators 
with r = 1, \p\ and k > 2 — r. It coincides with the statement of the table 1. 

Thus, we obtained that the modes of VF-constraint are positive or negative depending on 
where is the singularity of the potential. Certainly, it is rather natural and could be guessed 
from the very beginning. 

Let us say some words to justify the procedure we have proposed. We already stated in 
section 2.3 that string equation operator A (2.36) along with reduction condition imply the 
presence of constraint algebra which fixes unambigously (as perturbative series) the partition 
function. We obtained in this section all constraints which are necessary for this purpose. 

3 The BGWM versus the GKM 

After the presentation of generic theory of the GKM, we are now equipped to discussion of the 
BGWM as a particular example antipolynomial GKM. However, to begin with we still need 
to give a little more comments on the reasons, why the BGWM can indeed be identified as a 
GKM. In making this identification we shall be naturally lead to the introduction of a concept of 
"universal" BGWM, which in a certain sense is unifying such models for all the unitary groups. 

3.1 Ward identities for the BGWM and their GKM-like solution 

Integrability of matrix models is usually a corollary of huge covariance of the (matrix) inte- 
gral, which is used to define the full partition function. Since the action is of generic type, 
arbitrary change of integration variables results into some transformation of coupling constants, 
and invariance of the integral under such change implies restrictive constraints on the functional 
dependence of partition function on the coupling constants (external fields). Some of constraints 
are explicitly resolvable, but others form less trivial closed algebras. When these are isomorphic 
to (subalgebras of) some natural cartanian-type algebras (like Virasoro, W- or W-), solutions to 
the constraints are naturally r-functions of conventional (cartanian) integrable hierarchies (i.e. 
of (multicomponent) KP or Toda type). The set of constraints can be considered as invariant 
description of partition function, of which the original matrix model is nothing but particular 
integral representation. Other representations can differ by choices of integration contours and 
provide a kind of analytic continuation of the original function. See [1] for more details. 

According to this description, one should begin analysis of integrability structure, if any, 
of a given matrix model from identification of the adequate changes of integration variables 
and derivation of the corresponding constraints (Ward identities). The next step should be the 
choice of the coupling constants which brings these Ward identities to some standard form. The 
last step - identification of integrable structure - is yet not always possible, because the theory 
of non-cartanian hierarchies (i.e. those, not associated with the level k = 1 simply-laced Kac- 
Moody algebras) is not worked out in any detail. This, however, will not be an obstacle in our 
first example of the BGWM, which appears to suit into the standard Toda-lattice pattern. 

The BGWM partition function is defined by the integral (1.1): 

Z BGWM (J, J^) = tt~ ( [dU]e^ u+JU '\ (3.1) 

VN JNxN 

The "coupling constants" of the model are represented by the N x N matrix J (external matrix 
field). Among the admissible changes of integration variable U are left multiplications U — > VU 
by any unitary V, which leave the Haar measure [dU] invariant. The associated Ward identities 



read just 



2bgwm(J, J ] ) = Z BGWM (JV, VU^), 



(3.2) 



and together with their analogues, reflecting the right-multiplication invariance of [dU], they 
imply that is in fact a symmetric function of only N variables: eigenvalues rrii, i = 1 . . . N 
of the matrix M = JJ^ 10 : in particular, 

Z BGWM (J,J^) = Z N {M). (3.3) 

Dependence on these remaining variables is defined by the more involved Ward identities, which 
can not be resolved in such a simple way. It is most convenient to write them in the form of the 
matrix- valued identity [3, 4] 11 : 

d d 

— r • -^-Z B gwm(J, J^) = I ■ Z-qgwm{J, J^)- (3.4) 

0«/ t ' r OJtr 

This relation holds just because the derivatives at the l.h.s. produce a product U • U' = I 
under the integral. If Zjy(M) is now substituted instead of ^bgwm(^) J ) into (3.4), we get (see 
Appendix): 

ak M ak ZNiM) = I - ZMM) - (3 - 5) 

This is the equation that can look somewhat familiar from the theory of the GKM. Kontsevich 
integral (2.1), 

jr^jsf, L) = ±- I dXe ^ LX - m ^ x+v ^ (3.6) 

'n Jnxn 

is defined with the "flat" Hermitean Haar measure dX = Wa^dX^. This measure is invariant 
under conjugation X — > VXV^ with any unitary matrix V, thus 

F v (L)=F v (VLVi), (3.7) 

what implies that Ty is a symmetric function of eigenvalues l a , a = 1 . . . n of L only. 

Remaining less trivial Ward identities are associated with more general changes of integration 
variable. Of interest for us will be implication of particular transformation: X — ► X + XeX 
with some infinitesimal (X-independent) matrix e. Invariance of the integral for any e implies 
the following matrix-valued equation: 



d T d . d ( d \ 2 . ( d 



dL tI dL tI ' dL tI \dL tT J \dL tT 



Fy(N,L) = 0. (3.8) 



This equation becomes identical to (3.5), provided M = L, thus n = N, and 

M = n, V(x) = -. (3.9) 

x 

Eqs.(3.5) and (3.8) define Z^{M) and Ty{L) unambiguously, thus equivalence of the equations 
implies the identity (1.2), 

Zbgwm(J, J ] ) = Z N (M)/Z N (0), M = Jjt, 

1 r i ( 3 - 10 ) 

Z N (M) = F 1/X (N,L = M) = -L dXe Tr(MX-Nlo g X + ±)_ 

Vn Jnxn 



10 This conclusion would not be true for SU(N) integral. In that case V € SU(N), thus DetV = 1 and Zn is a 
function of M = Jj\ and also of DetJ, DetJ^. 

11 Subscript "tr" here and below denotes transposed matrices, / stands for the unit matrix. 



The last relation can be also rewritten in terms of Y = 1/X: 

Zn (M) = -1- / dYe ^±-Nlo g Y + Y) (311) 
VN JNxN 

where we used the transformation law for the measure, 



1 dY 
Y ~ (detF) 

and thus 12 



This property in fact imples the symmetry of the integral (3.10) under the change of integration 
variable X <-> 1/X. Together with left- and right- invariance of (dX), 

(d(GX)} = (d(XG)} = (dX) for any G, (3.14) 

it can be used to restore the symmetry between J and jt (which is not really explicit at the 
r.h.s. of (3.10)): 

Z BGWM (J, J^) = ^rl dXe^ JX - N ^x+-^) (3.15) 

VN JNxN 

Invariance under X <-> l/X now implies the required identity 

Zbgwm(J, = Z B gwm(J\ J)- (3.16) 

3.2 Virasoro constraints and integrability in the character phase of the BGWM 

We are now in prepared to continue discussion of the BGWM, considering it as a particular 
example of the GKM with potential V(x) = 1/x. 

3.2.1 BGWM in the character phase 

In this limit partition function, 

2BGWM(^,^) = _ J nW 

Z N {M - 0) ^ 

z n (m = jjt) = r 1/x {N, jjt} = z+{t+)\ tt=1 _ w)k , 

is expandable in a series in positive powers of J, jt-fields. It can be also considered as a generating 
functional of the symmetric unitary matrix integral with Haar measure: 

^ + ) = i + E ( E c »M «i''' k £H) > ( 3 - 18 ) 

j£>i \i<k^.<k M (h + ... + k M y.) 

where the coefficients CN{k a } are defined as: 



To emphasize the relation of (dX) to Haar measure, note that for X = , we get the standard expressions 
for Haar measure [20], (D ** )N = Det ^f +g2) . 



±f[dU]UvUt = 6«6 j3 c N {iy, 

i- J [duPijUuUt.fi = (s^S^ + S^SfyS^) c N (l, 1) + (<VVV>/, + W«) ^(2); 

and so on. 

(3-19) 

The series (3.18) is naturally graduated by the number of U — U^ pairs in the correlator, which 
is equal to K{k a } = J2k a - Coefficients CN{k a } are defined through the recurrence relations, 
which are nothing but implications of the Ward identities for the BGWM. In order to derive 
them in explicit form we need to substitute Z^(t + ) into eq.(3.5): 



3 -M 1 £ r Z+(t+)=I-Z+{t+) 



dM tT dM tI 



t+ = \TrM k 



(3.20) 



and use (2.40) to rewrite it in terms of t-variables (see [1] for comments on this type of deriva- 
tions). The result reads: 



]TM S (w£i 2) (t)-5 s , )z+(t) = o. 



s>0 



(3.21) 



Since Ws + ' 2 \t) = Cs a> (N,t) are just "discrete- Virasoro" operators, 



s-1 



4 d) (^) = A^ + E^7^ + £ 



d 2 



dt s k>0 dt k+s k=l dtkdt s ^ k 



, s>+l, 



we obtain a set of discrete- Virasoro constraints for Z^(t + ): 

cW(N,t)Z+(t) = 6, tl Z+(t), s>+l. 



(3.22) 



(3.23) 



These provide the required set of recursive relations for cn {k a }, which can be used to derive the 
somewhat non-trivial explicit examples. If 



K{k a }-\ 

C N {k a } = C N {k a } [] (N 2 -l 2 )-\ K{k a } = Y / ka, 
1=0 

then cjv(&i • • • fejw) are polynomials of degree M in N, e.g. [21] (see also (2.59)-(2.60)): 

c N (l)=N; 
c N (2) = -N, cjv(1, 1) = N 2 ; 
c N (3)=4N, c N (l,2) = SN 2 , c N (l, 1, 1) = N{N 2 - 2); 
cjv(4) = SON, c^(l,3) = +8(2iV 2 -3), c N (2, 2) = +3(iV 2 + 6), 
cjv(1, 1, 2) = SN(N 2 - 4), c N (l, 1, 1, 1) = iV 4 - 8iV 2 + 6; 



(3.24) 



(3.25) 



Alternatively, the same quantities can be directly derived from the GKM-representation of the 
BGWM. 

The fact that poles occur in CN{k a } at all N < K{k a } is referred to as the De Wit-t'Hooft 
anomaly [22]. These describe the singularities of Z~^(t) at generic values of t-variables. This 



function is a kind of a universal object, describing all the BGWM models (for all unitary groups) 
at once. Reduction to particular group implies that time-variables are substituted by the N- 
dependent quantities tjj" = |Tr( JJ^) fc , so that only N of them remain independent. Restriction 
to this hypersurface in the infinite-dimensional space of time-variables spoils some nice proper- 
ties of Z~x(t), but instead on this hypersurface the De Wit-t'Hooft poles are canceled between 
different terms with the same K{k a }, in accordance with the finiteness of the unitary group 
integral (3.1). 



3.2.2 Universal BGWM in character phase 

In order to establish connection with the pertinent GKM it is necessary to find a representation 
of the BGWM partition function in terms of an integral, where the size of the matrix is a new 
parameter, independent of N. The answer to this question depends, of course, on the way 
how the dependence of N is separated from that of all other variables. We choose eq.(3.18) as 
the definition of the universal BGWM in the character phase (i.e. t + 's and N are considered 
to be independent variables). Invariant definition of this quantity is provided by the Virasoro 
constraints (3.23), and our goal now is to find a solution to this equation in the form of an n x n 
Hermitean matrix integral with any n. This is a simple exercise if we return back to eq.(3.8) 
and make the substitution V(x) = 1/x, but do not identify N with n and N. Then we have: 

(note that according to eq.(2.1) T is defined as an n x n integral). Let us now consider 
Fi/x{N ,L) = Z{t + \N) as a function of t\ = ^tiL k , k > 1. There is no explicit n-dependence 
in Z\ It appears, however, in eq.(3.8). Indeed, the first operator at the l.h.s., when acting on 
Z, turns into J2s>o L s £^+i(n,t + ), where n appears explicitly. Moreover, the contribution from 
the second operator turns into (n — J\f) J2s>o L s df +1 ancj - ^ oes not cance l this n-dependence, 
instead it changes C^ d \n,t) into C^ d \2n —J\f,t) and the resulting constraints are 

£ { s d) (2n - M,t)Z + {t + \M) = 5 s ^Z+{t + \N). s > +1. (3.27) 
We see now that (3.23) can be reproduced, if N = 2n — J\f, i.e. N = 2n — N and 



Z nW , uju= ^M 2n ~ = I dXe^ LX+ ( N - 2n ^ x+ ^. (3.28) 

Of course, for n = N we obtain the old formula. Eq.(3.28) can be also rewritten in terms of 
Y = 1/X, thus getting a GKM with non-standard coupling of the L-field: 

Z+(t)\ HLk = ^l dYeW 1 -****^. (3.29) 

In this form it can be interpreted as a generating function for correlators in the Penner model [7] . 
This expression also reminds the theory, discussed in ref.[8], though there are some important 
differences in our interpretation of this integral. 

Since Z^(t) satisfies n-independent Virasoro constraints (3.23), it does not actually depend 
on n. 13 This n- independence is crucial for interpretation of Z^(t) as a KP r- function. 

13 This is despite explicit n-dependence of the action in (3.28)! One should first express everything in terms of 
t- variables and then observe the elimination of n-dependence. The issue of n-independence is somewhat obscure 
in ref.[8]. It is not very clear whether it can be preserved for generic potential V(Y) / Y. 



3.2.3 Zt T as a r-function 



We proceed now to determinant representation of the integrals (3.28), (3.29). We derive explicit 
representation in terms of Bessel functions and demonstrate the n-independence explicitly. 

Application of the usual GKM technique to the case of (3.28) gives: 



Jnxn J _■ 



A(Z) 

-det^ / n^-^ 1 /, = det *»fo) 



where the would be "basis vectors" are: 



tf o (0 = y ^ x «-l e ^-^log^+V^ = J dy y -l-a e y+Nlo g y+l/y = ^ ^ j dy y -l 

k=0 



(3.31) 

where normalization of ^ a (0 can be left unspecified at the moment, I a {z) are modified Bessel 
functions and (3.31) is correct for all real values of M (in particular, there is an evident symmetry 

Ia = I-a). 

To discuss the dependence of n we make explicitly the steps which were explained in section 
2.2.1 for the general case. 

The basis vectors correspond to character phase, therefore, their asymptotics (at small U) 
are ^/ a (0 ~ 1 + 0(1). This means that one should do the transformation / — > A -1 and consider 
linear combinations of a 's not changing the determinant (3.30) as it explained in section 2.2.1 
to adjust these linear combinations such that new vectors * a (0 will have the asymptotics l a ~ l . 
Let us also normalize them to have the unit coefficient in the first term. 

So far the problem of the normalization was left beyond our consideration. In fact the Ward 
identities discussed in the previous subsection knew nothing about normalization. Therefore, 
we are free to fix it arbitrary. In accordance with (1.1) we prefer to use the requirement that 
Z = 1 at all Zj = (or all = 0). This is achieved by adjusting the coefficients in front of the 
leading terms in the small-Z expansion of * a (0 to unity. 

With these requirements in mind one can derive the following formula for 9 a . 

9(1)- V lk (3 32) 

" JZ-i r ( a — -A/* + k + l) T(k - a + 2) ' (3 ' 32) 

In order to allow comparison with the r-function theory let us substitute I — > 1/A. It remains 
to rearrange the indexes, a — > n — a + 1, and obtain: 



where 



_ dety> a (A b ) 

TtL*=£.Ar fc - A (A) ' { ' 



, ,^ T ,u T „ ^ xav- T((2n-AA)-2a + 2) \~ k xaV ^ fc 

(3.34) 



and matrix p a i~ can be continued to the non-positive values of k by definition p a j, = 0, whenever 
k < (this definition easily follows from the manifest form of p a k (3.34) due to the poles of 
T-functions in the denominator). 

Thus, we can see that to eliminate any n-dependence of the partition function (3.33) we can 
just choose M to be 2n — N, where N is a free parameter. Indeed, this parameter N is nothing 
but the size of matrix in the BGWM and, simultaneously, it appears to be the zero time (times 
2), as can be understood from (2.9). 

Now let us rewrite the sum (3.34) in terms of modified Bessel functions and in integral form 
like (3.31): 



r(iV-2a + 2) 
r(iV-2a + 2) 



A 



a-l 



dyy 



2a-N-2 e l/(Xy)+y = 



2 2-a 




N-l 



(3.35) 



lN/2-ai^). 



One can display manifestly in the expression (3.34) the dependence of zero time, in accordance 
with (2.9): 



Z+(t)=T N/2 (t) 



_ det ab X N / 2 ^ N/2 (0,X) _ det ab X N /^ a _ N/2 (X) 
^=1 Ei XT" ~ A(A) = A(A) • (S - 6b) 



Therefore, we obtained that in the character phase the partition function of the BGWM is the 
r-function of KP hierarchy, corresponding to the element of the Grassmannian given by the basis 
vectors (3.34) with N/2 being zero time. These basis vectors give rise to the reduction of the 
following type. Let us consider infinitesimal additive variation of variable y in (3.35). It induces 
the recurrence relation for the integrals of the type of those in (3.35). It connects integrals with 
even and odd degrees of y in the integrand. Applying the procedure successively three times we 
finally obtain the relation: 

* a = A ^" 1 + 2(a - l)(a - 2/ - 1 " (2a-3)(2a-4)^(2a-5) V ' a - 2 - (3 ' 37) 

This is a counterpart of the reduction condition (2.11) in the character phase which is, in fact, 
a sort of Toda chain-like reduction (see sect. 2 of the paper [10]). 

From the expressions (3.33)-(3.34) one can see that, indeed, this partition function is nor- 
malized to be unity when all times are equal to zero, or all Aj — > 00. 



3.2.4 Character representation 

Let us now demonstrate how it is possible to use these explicit expressions to obtain the coef- 
ficients Cjv{fc a }. The expansion coefficients Cjv{m a } of a generic GKM in terms of the unitary 
group characters [2] can be written down for the particular case of (3.28) and comparison of the 
formula 

00 

z nW = E C N {m a } X mi ,m,2,... (*) (3.38) 

mi>m2>...>m„>0 

with (3.18) provides an expression for CN{m} in terms of CV{m}. 

The first step is to work out (3.38) in explicit form. Using the standard formulas from [2] 
(see subsection 5.4 of that paper), one obtains: 



mi>m2>...>m„>0 vr / 



^ det p a ,m b +a-b+l 



00 T(N — 2a + 2) 



(3.39) 



where 



det 1 \-m b +b-l 

This is the first Weyl formula for (primitive) characters, A a are the inverse eigenvalues (in 
contrast to the notation of the paper [2]) of a given unitary matrix, and {m a } are the lengths of 
the rows of corresponding Young table (therefore, they should be ordered as in formula (3.39)). 

Now, one can see immediately from (3.39) that, if some m r+ i = (and, therefore, so all 
m q , q > r), then P{rrii} will be determinant of block matrix: 



-P{m a } = det 



A 



D 



C 



det A = det Pa,m b +a-b+l 
rxr ab 



a,b<r 



(3.41) 



Thus, we obtain that the final answer for the expansion (3.39) does not depend on n. The only 
restriction to n is that it should be larger than the number of non-zero m a 's, i.e. there should 
exist such p < n that the representation with these m a 's can be embedded into the group U(p). 

For integer N formula (3.39) can be transformed into expression 



mi>m2>...m r >0 
1 

ril<a<6<Af (l + m b -T b ) 



(mi !m2! • • • m r !) 2 



n 1 



Ka<b 



m a + b - a 



(3.42) 



(*) = E C N {m a }x mi , m2 ,..., mr {t) 

mi>m2>...m r >0 



which coincides with that in the paper [21]. The only difference is that author of [21] considered 
./V equal to n, and N, therefore, was not just parameter, but the size of matrix. It explains why 
his answer has no De Wit-t'Hooft poles, while our universal function, which can be determined 
from the coefficients of expansion in the characters (see below) does have. Indeed, the term in 
the first brackets does not depend on N at all, while the term in the second brackets is nothing 
but inverse of dimension c/jv(mi,m2, ... ,m r of the representation of the group U(N), which is 
given by the numbers mi, m2, . . . , m r . The simplest representations have dimensions: 



(3.43) 



tt (N + m a — a)\ f-r 

d N (m 1 ,m 2 ,...,m r ) = 1_1 + r - a)!(iV - a)! X 11 K> - m b + a - 6). 

This expression, and, therefore, the second factor in the formula (3.42) can be trivially continued 
to non-integer N replacing factorials by proper T-functions. Then, one can trivially see the 
presence of poles at all integer values of N = 0,1, . . . ,r — 1. 

Thus, in the universal function, when considering N as a free parameter with no the reduction 
to the hypersurface in the time space corresponding to finite n = N, there are De Wit-t'Hooft 
poles. 

Now let us say some words on how one can effectively obtain from the expression (3.39) 
the coefficients C7v{A; a } of the time-expansion (3.18) of the universal function Z^. For further 
convenience let us rewrite the coefficients C7v{m a } in (3.42), using (3.43), in the form 

- n r( W+ J- w a riK,'! +r - a )! x n<-. --».+.-»), 

where again only the first product depends on N. The first coefficients are equal to 
CW(1) = jf', 

Cn{2) = 2 jV(jv+l)' Cn(1, 1) = 2 N(N-1) ' 

CV(3) = 6N (N+l)(N+2) ) CW(2, 1) = 3AT(Ar + i)(jV_i) > CW(1, 1, 1) = 6AT(JV-l)(iV-2) ! 

{CV(4) = 247V(7V+l)(Af+2)(7V+3) ' C< A r ( 3 ' *) = 2iV(iV+l)(JV+2)(iV-l) ' ( 3 ' 45 ) 
Cjv(2, 2) = 12 iv a (JV+l)(iV-l)' C^v(2, 1,1) = 8iV(iV+l)(Ar-l)(iV-2)' 
Cjv(1, 1, 1, 1) = 24iV(iV-l)(JV-2)(iV-3)(Ar-4) 5 

The next step is to rewrite the characters in terms of times. To do this, let us introduce Schur 
polynomials Pk(t) defined by the expansion: 

exp{^ fc x fe } = £P fe (t)z fc , k<0, (3.46) 
the polynomials with negative indices being put zero by definition. Some first polynomials are 



Po(t) 
Pi(t) 

W) 
W) 

p 4 (t) 



1, 

h + ^ 



2 ' 

t 3 + ht 2 + t i , 

2 



(3.47) 



U + hh + | + ^ + ^ 



4 
1 

24' 



Now we can use the second Weyl formula for (primitive) characters 



(t) = det P ma . a+b (t) (3.48) 

ab 

to rewrite characters through time variables. Now we can calculate the coefficients C]\r{k a } 
immediately using the formulas (3.38), (3.44), (3.46) and (3.48). Say, using some first characters 

X2{t) = t 2 + i, Xii = t - *a; 

^3 ^3 ^3 

X3(*) = *3 + *i*2 + X2i(t) = -j -t 3 , Xm = h - ht 2 + -5-; 
' X4 (t) = t4 + tlt3 + 4 + ^ + |, X3l(t ) = ^-t 4 -4 + |, (3.49) 

< X22 ( t ) = 4 - ht 3 + ^, X211 (t)=t 4 -4_^ + 4, 

Xim(*) = *i*3 - *4 - ^ + t + 
one trivially reproduces the expressions (3.25). 



3.3 Virasoro constraints and integrability in the Kontsevich phase of the 
BGWM 

3.3.1 The BGWM in the Kontsevich phase 

We need now to consider eqs. (2.3-2.6) with L = M and n = N for the particular case of 
V(X) = y and M = N. In such situation Xq = M~ x l 2 and p = —2. The C-factor, appearing in 
the definition of partition function Z~ (i~(~ 2 )), 



^BGWM(J,J t )= Z ^ = J J ) ) , 



Z N (M)=F 1/X (N,M) = C 1/x (N,M)Z-(t^) 
is given by (2.6): 



(3.50) 



-2fc+l 2fc-l Lim 



Cl/X ~ (DetM)- JV / 2 Det 1 / 2 (M 1 / 2 ®M + M® M l l 2 ) ' 

so that 

^r ( ~ 2) )L-(- 2 )_ x TrM - fc+1/2 = ft V^^P^ 2E » my2 ^W. (3.52) 

Ward identities for Z~ result from the substitution of Z^{M) into (2.39) 14 . Though general 
procedure was already written in section 2.4.3, this is the only case which was worked out for 
antipolynomial model [6]. The result is the set of Virasoro constraints imposed on the partition 
function: 

£ s Z-(t- ( -V) = 0, (3.53) 



14 This substitution was already performed in ref. [6]. Since the relation between the GKM and the BGWM was 
not discussed in that paper, expression (3.52) for the C-factor was just guessed and postulated, without reference 
to generic prescription of the GKM. 



where 



1^.5 1 d 

(3-54) 

1 V- , 9 9 

fe<0 

Let us point out that, in accordance with the general rule (2.61), the shifted time is t-±, and 
the first constraint is Co. As it was already explained in section 2.3.2, this constraint is the 
string equation in P = — p = 2 case. Along with the statement that the partition function is a 
r-function of KdV hierarchy it is sufficient to fix the partition function, as well as to reproduce 
all the tower of Virasoro constraints (3.53). 

3.3.2 The Universal BGWM in the Kontsevich phase 

Now let us emphasize that the Virasoro algebra (3.54) does not contain N dependence at all. 
It is, certainly, the trivial consequence of the fact discussed in section 2.2.3, where we have 
demonstrated that the partition function of the BGWM is a (r)-function of times, which does 
not depend on N. It means that we do need no additional parameter to define properly universal 

function in the Kontsevich phase. Indeed, in this case the function Z~(t~^) _ (2) is just needed 

t k 

universal function which does not depend on the size of matrix. This is why we omitted the 
subscription N in the formula (3.52). 

Still some additional parameter can be introduced into the partition function, namely, zero 
time, how it was explained in section 2.2.3. This parameter is completely analogous to the 
parameter N in the character phase, but has nothing to do with the BGWM, as the BGWM 
corresponds to zero value of this parameter (see (3.52)). 

It stresses the difference between structure of the universal functions in the different phases. 
Generally, in the Kontsevich phase the partition function has more regular and universal behavior 
which reduces the number of essential parameters, unlike the character phase when each change 
in the potential requires a special treatment, therefore, having many parameters which govern 
the behavior of the universal function. It explains why the integrable treatment which gives, in 
a sense, a universal description, is hardly applicable to the character case. 

3.3.3 Z~ as a r-function 

As we already discussed in subsection 2.2.3 the partition function in general antipolynomial 
Kontsevich case is the r-function of KP hierarchy specified by the reduction condition (2.26) 
and string equation. It can be described as well by its explicit determinant representation (2.9) 
with the basis vectors determined in (2.24). 

In the concrete case of BGWM, one should put p = — 2 and get KdV r-function with basis 
vectors having the manifest expression (see (2.24)-(3.31)): 



^a-iv(A) = 2^Fe VAA J a _ A r_ 1 (2A) ^ A ^——^———^ + ... 

(3.55) 



15 In the reference [6] the signs of indices of all Virasoro generators and times were chosed opposite, and the 
numeration of times is a bit different. 



Parameter N here is just zero time. The dots in the asymptotics expansion stands instead of 
exponentially small at large A terms. 

This manifest expression for the basis vectors allow us, among other, to prove explicitely the 
independence of the r- function of even times, which was noted in [6]. Indeed, the reduction con- 
dition to KdV (2.26) already implies that the r-function can include any even-time dependence 
only as an exponential, linear in these times (see (2.12)). Then, we can use the trick proposed in 
[9]. That is, we can choose only two non- vanishing Miwa variables Ai = A and A2 = —A. Then 
all the odd times to be zero, and independence of the full r-function of even times is equivalent 
to the fact that r = 1 in this particular case. It is trivial to check that the r-function is equal to 



-2ttA 



{v / A/ jv (A)}{v^A/jv-i(-A)} + {v /z A/^v(-A)}{Va/ J v_i(A)} 



(3.56) 



in the point under consideration. The most subtle point is what we should understand by the 
function I c {—z) which is ambigously continuable to the negative values of argument. Indeed, 
how it is usually understood in GKM approach, we consider this object as continuation of formal 
power series. This means that we consider instead of I c {— z) Macdonalds function ^K c (z), where 
the normalization is chosen to have the same asymptotics as in (3.55). The only difference of 
these two functions as formal series is in exponentially small terms which are just governed by 
Stock's phenomenon and are not taken into account in any integable treatment 16 . Now we can 
use the identity between functions I c {z) and K c (z) [23] 



I c+1 (z)K c (z) + I c (z)K c+1 (z) = - z (3.57) 

to get finally that (3.56) is equal to unity. Thus, we have proved that the r-function of p = — 2 
model does not depend on even times at all. The proof for the general case can be done mostly 
like that proposed in the paper [9]. 

Now we would like to make some comments on the literature concerning this r-function. To 
begin with, let us note that any product of two KdV r-functions taken at sucsessive values of 
zero time (say, and 1), is an MKdV r-function (see, for example, [24, 25, 5, 26]). The MKdV r- 
function obtained as a product of two p = —2 GKM r-functions 17 can be specified by constraint 
algebra (or, sufficiently, by its lowest constraint, i.e. by the string equation) (3.53)-(3.54). As 
the constraint algebra does not contains any A-dependence, we can write 



C q T s = 0, q > 0, 

(3.58) 

C t s = 0, 

where s = 0, 1 labels r-functions with two sucssessive values of zero time. 

In the papers [5, 26] it was demonstrated that the general partition function of unitary matrix 
model [27] (not in the external field, but with arbitrary potential) in the double scaling limit is 
MKdV r-function which is the product of two KdV r-functions satisfying the string equation 18 



C q T s = 0, q > 0, 

(3.59) 

16 These terms can effect, in particular, the possibility to present the determinant of basis vectors as a function 
times. 

17 One of them can be chosen to be BGWM partition function, but the other one can not be, as there is no way 
to introduce zero time to unitary matrix integral. This point caused a problem in ref. [6], where an attempt to 
construct the BGWM unitary matrix integral for MKdV r-function was made. 

18 Their operator £0 does not contain the constant 1/16, therefore, our definition of the parameter ^ is a bit 
different. 



all other constraints being just Virasoro algebra (3.54). fj, here is just a free parameter which 
is not specified at all. Comparing this expression with (3.58), one can see that p = —2 GKM 
gives an explicit example of the double scaled partition function of the unitary matrix model with 
\i = 0. To construct the general case, one just needs the KdV r-function which satisfied the 
deformed string equation (3.59). 



3.3.4 String equation 



In section 2.3 we derived the string equation in some non- manifest way. The different approach 
to this is to use the direct derivation described in the very beginning of subsection 2.3.1. This 
derivation is less useful and a bit tedious but it can also be done. We demonstrate this for the 
simplest case of p = —2. 

Indeed, let us consider the derivative of the BGWM partition function (cf. with (2.27)): 



Tr 



A 



d 



dA 



tr 



logZ-(H- 2 )). 



(3.60) 



It can be written as 



ETr 

p<0 



A 



dt^ 
dA 



p+i 



tr 



dt 



d 

2p+l 



1 ' 



logZ-(t-(- 2 )) = -]TTr 
J2 (2p+l)t 2p+17 J— log Z-(t 



dt 



d_ 

2p+l 



log Z-(t-^) 



p<0 



2p+l 



(3.61) 

On the other hand, the derivative in (3.60) can be manifestly calculated. To do this, let us 
represent the BGWM partition function in the determinant form (2.9), (3.55) 



2k + l - 



det ip a (\ b ) 
A(A) 



with 



(3.62) 



^ (A) = e~ 2X V\ J dxx~ a e X2x+l ' x = e~ 2X F a (\). 



(3.63) 



Now we need two facts which are true in the general antipolynomial case, i.e. our proof can 
be immediately generalized. 



Let us note that the action of the operator A 2 ^- transforms F a (X) into F a+ \ (A)+ J2 FbW 



Kb<c 



(we used the reduction condition (2.26)), i.e. there are linear combinations F a (X) of F a (X) such 
that F a+ i(A) = X 2 -^jF a (X) = (a 2 ^-^ ^i(^)- The second important property which is crucial 
for our derivation is the asymptotics behaviour of F a (X) at large A (see (3.55)): 



Fi(X) - V 4 

Now we are ready to calculate the derivative (3.60): 



- A" ( I - N ^ N . ^ + oil A; ) for any a > 1. 



(3.64) 



Tr 







OA 



tr 



log Z-(r { -V) = ~2J2 A a - E A a ^-logA(A)+ E A a ^-logdetF a (A) = 



= "2 E A* + + 2 S A ^ logdetF a (A). 

a a ci 

(3.65) 

To calculate the last derivative we use the trick from the paper [9], where it was demonstrated 
that (the derivative with respect to the first time in [9] is to be replaced by that with respect to 
t-i here) 



dt-! 



logZ~(H- 2 )) = res A 



^i(A) 



(3.66) 



L n j= i(A-A,) 
" iE A a^|logdetF a (A)[l + 0(l/A)]^ 

if the conditions (3.64) are satisfied along with the determinant representation (2.9) with the 
entries (A 2 ^) a Fi(A). Therefore, using (3.64) we obtain 



d 



^ log^r ( - 2) ) = ± - -E^logdetF a( A). (3.67) 

Now, collecting together (3.61), (3.65) and (3.67) one can finally obtain the string equation: 



+ 2^- + -\z-(t-^) = 0. 



(3.68) 



This is evidently Lo-constraint which coincides with the formula (3.54). 

It is interesting to note that in the course of the derivation there arised the manifest de- 
pendence of n, but it disappeared from the final result. This reflects the rather complicated 
n-structure of antipolynomial Kontsevich phase. 



4 Itzykson-Zuber integral as a 2-matrix model 
4.1 From the BGWM to IZ integral 

Since partition function Zn(J,J^) = Zn(M) is a generating functional for any correlators of 
{/-matrices with Haar measure, it can be used in particular to represent the IZ integral in terms 
of Hermitean matrix model: 

Z IZ ($,$|J, Jt) = J- / [dU]e ^U + JU%TrZU^ == . Z BGWM (J,$-^-d> + Jt) : = 

d _ If e Tr ^dX _ Tr(j*^-+jjt)x _ 



: Z N { JJ^ + J^jpr®) : = ^- I 
oJtr V N J]\ 



V N JnxN (BetX) N 

(4.1) 

"Normal ordering" sign at the r.h.s. implies that all the derivatives d/dJ should be placed to 
the left of all J's. Pulling back all these derivatives to the right provides the following relation: 



(4.2) 



If we now put J = J< = the answer arises for the IZ integral per se: 19 

*) - ^ / [cH7]e™^ = / ? #^e^>o (4.3) 
Vn JnxN JnxN (DetA) JV 

We can now make a change of matrix variables 



X -> H = — . (4.4) 



Invariance of the measure implies that 

(DetX) N (Dettf )^ ' 

It remains to denote 

T+ = ^Tr<§ fe , fe > 0, 
= Tr I = N 



(4.5) 



(4.6) 



fc=0 

(compare with our usual definition of the "+"-time variables) in order to get: 

Vat JAfxTV Vat JAfxAf 

. (4-7) 

This represents IZ integral in the form of the GKM with generic potential V(x) = J2k^k x k - 
The first representation is identical to expression from [8] for the (non-full) partition function 
of the c = 1 model. 

One can also work out similar representations at large, rather than small eigenvalues of 
It is enough to note, that 

Det($® J-J®ff) " ' ( j 



where 



1 I <l> ". / > IS. 

(4.9) 



T~ = jTv^- k , k>0, 



= -iVTrlog$ 

k=0 



Thus (4.7) can be also represented as 



Vat JAfxAf Vat jatxAt 

(4.10) 

These expressions are, however, asymmetric in $ and <!>. Symmetry is restored in another 
representation of IZ integral - in terms of conventional 2-matrix model. 

19 Another version of this formula (for the derivation, see Appendix): 

1 f c Tr(i-iVlogX) 



Vn Jnxn Det(/ ® / - $ ® $A) 



4.2 Relation to 2-matrix model: direct proof 



Now we are going to state the connection (1.8) between IZ integral and two-matrix model by 
direct calculation. We need only the formula for the Cauchy determinant 

v Xi-yj Ui,j{xi-yj) 

and the integration is chosen such that 

f j f{h)f{h)dhdh = J f(h)dh J f(h)dh. (4.12) 
Now one can write the chain of identities: 



™ I I dHdHe™»e^>-^ HTk+T ^ = 

V N J J NxN 

M r r e TrHR dHdH 

Vt J JnxN BeUH ipetlff i I - / i 5) ~ 

, ^ ,^ (4-13) 

1 f/n. JVe hA A(h)A(h) (ui) 



y y n 



V!./ ■/ 11 H,.,i/<, o ; ) 
' f TT dhidhie hi ~ h > det — ^ det - 1 - = - 1 ,- s det , 



AnA(0)A(0) i i Y ij /»i - ^ ij ^ - ^ A(0) A(0) 

and we used at the last stage the convention (4.12) which, after manifest writing the determinant 
as sum over permutations, leads to N\ equal terms of each type (of each possible permutation) 
and transforms to N\ equal determinants. 

Now let us say some words on the integration contour in these expressions. In fact, it is 
different in the cases of potentials, which are polynomial and antipolynomial in H. Indeed, in 
the first case we use the expansion of logarithm in the domain where hi/4>j < 1 for any i and 
j (see the first equality in the chain (4.13)). It means that in the course of integration (last 
equality in (4.13)) all 0j's lie out of the integration contour, or, equivalently, the integration 
goes around oo. Unlike this case, antipolynomial potential implies <j>i/hj < 1 in the first equality 
in (4.13), i.e. integration contour is closed around zero, all <fo's being inside it. Under these 
rules, the integration is to be determined as usual one complex variable integration over closed 
contour: 



/ 



-dh = e^ h (4.14) 



h 

and so for the integration over h. It completes the proof. 

In the next subsection we discuss some properties of the IZ integral which are consequences 
of this its representation in terms of two-matrix model. 



4.3 IZ integral as a r-function 

The real point is that we know many different facts about two-matrix model, and, among others, 
that its partition function is a r-function of lattice Toda hierarchy. Therefore, we can assert 
that the IZ integral is a r-function in Miwa parameterization (1.9). It means that we can fix 
the point of the Grassmannian which corresponds to this r-function and obtain some useful 
representation just using the theory of integrable hierarchies. 



The point of the Grasssmannian can be described in two different ways. The first one is to 
express it through exponential of a bilinear of fermions. It was done in [2]. Another possibility is 
less manifest, but is directly connected with the formalism of two- matrix model. Namely, using 
the formalism of orthogonal polynomials, one can rewrite (1.8) in the determinant form [28]: 



(4.13) = / / n^^e' l<S<+E ^ o(Tfc± ' ir+ ^ ± ^ fc) A(/ i )A(^) = det Hi j(T, T) (4 .15) 



Nl J J 11 V ' V 1 NxN 

with 



Hi 3 \ 



,(T,T) = / M^ e i&fE^+W (4.16) 
and properties 



(4.17) 



In fact, the determinant representation (4.15) along with the properties (4.17) give rise to a 
generic r-function of Toda lattice hierarchy, and the formula (4.16) specifies it. 

There are different possible applications of the integrable structure of the IZ integral. Say, 
having the formula (4.15) one can trivially find the fermionic representation for the element of 
the Grassmannian [10] and, then, by the Miwa transformation of times restore the IZ integral 
(instead of the manifest calculation of the previous subsection). Certainly, it can be also done 
reversly. But now, as an example, we are interested in only one possible implication of the 
integrable properties. Namely, we can find out the expansion of the IZ integral at small times 
like it was done in subsection 3.2.4 in the BGWM. Moreover, it can be done equally for both the 
polynomial and antipolynomial potentials in (1.8). This means that one can equally construct 
large- and small- (0, 4>) expansion. 

For the definiteness, let us consider the antipolynomial case. Then, the integration rules of 
the previous subsection implies that 

/ x a y l3 e xy dxdy = } - ■ (4.18) 
J r(— a) 

Then, using (4.15), one can obtain the determinant representation in terms of (Schur poly- 
nomials of) times 20 : 

■ let e**' = det £ — / — P t (^(^, a (4.19) 



AU)AU) i,j run ^ T(N + i - n + 1 



where the sum over i can be extended to all integer, the real range of summation being deter- 
mined by poles of the r-function and by vanishing Schur polynomials with negative subscriptions. 
For example, the first term of expansion of antipolynomial case looks is 

20 Certainly, one can trivially write down an analogous expansion also in the polynomial case. 

21 Note that in generic case (i.e. when formula (4.16)) is replaced by a general solution Hij(T,T) = 

J dhdhfi(h,h)e'^' k i k k with an arbitrary measure function n(h,h)) r-function has the analogous 

expansion tn(T,T) — detij m Pi- i (T)Ti m P rn -j(T) with an arbitrary matrix T t j = J dhdhh l h j fi(h,h) which 
describes the adjoint action of the element of the Grassmannian on fermionic modes. 



detV — — -Pi{T)P m+i JT) = , T , 1 (1 + . . .) • (4.20) 

Generally speaking, the calculation of the expansion is more complicated in comparence with 
the BGWM. The underlying reason is that now the number of Miwa variables A can not be 
made independent of the size of matrices in IZ integral. 



5 Conclusion 

This paper defines a framework in which both the phase- and integrable structures of matrix 
models can be discussed simultaneously. The crucial role is played by the concept of universal 
partition functions (UPF). Each UPF describes a set of matrix integrals, differing by the size n 
of the matrices which are integrated over and serve as external fields. UPF is defined to be an 
n-independent functional, depending on the form of the action (i.e. on the choice of the matrix 
model), but not on the matrix size. 

This definition of UPF depends on the choice of the variables ( "coupling constants" ) , which 
are considered to be properly reducible under the certain embedding of the smaller n\ x n\ 
matrices into the bigger x n2 ones. For the simplest embedding, when the smaller matrix 
is supposed to stand at the left upper corner of the bigger one (which is implicitly accepted in 
this paper), the example of the properly reducible variables is provided by traces of (powers of) 
matrices. 22 For eigenvalue matrix models this is also a complete set of variables. This means 
that whenever partition function of a matrix model is expressed in terms of the infinitely many 
"time-variables" t^ {p) = ±tr nxn M ±fc / p , which are treated as independent variables (despite only 
n of these are actually independent), and there is no explicit n-dependence in the shape of this 
function Z{t}, it is properly reducible and can be considered as UPF. For numerous interesting 
matrix models such UPF can be further identified with the r-functions of integrable hierarchies 
(if matrix model is of eigenvalue type, these are usually multicomponent Toda,- i.e. "Cartanian" 
- hierarchies). 

The best way to specify UPF is through the system of equations, of which the Virasoro- and 
W-constraints are the simplest examples. Original matrix models are then particular solutions of 
these equations, represented in the form of n x n matrix integrals, which describe the restriction 
of the UPF on peculiar n-dimensional hypersurface in the t-space. 

Construction of the set of the n x n integrals, possessing the same UPF is a matter of art. 
Moreover, given A x A integral can be associated with different UPF - depending on the par- 
ticular choice of the time- variables (even if the abovementioned embedding of matrices is fixed). 
We illustrated this phenomenon in the main body of the paper by consideration of the same 
Kontsevich integrals in two different asymptotics (character and Kontsevich "phases"), where 
the adequate time variables are either i£ or t k ^ p \ The corresponding universal partition func- 
tions are not the same! (For example, Z^(t + ) from section 3.2 depends on A, while Z~(t~^) 
from section 3.3 does not.) All these notions and features are most transparent in the study of 
Generalized Kontsevich model, which is the basic one in the entire field of the eigenvalue models 
and was in the center of discussion in this paper. 

Above scheme can help to reveal the relation of various physically relevant U(N) lattice 
Yang- Mills models to integrability theory. Partition functions of these theories are often strongly 
dependent on the parameter A, and this dependence is physically relevant and, if present, does 
not need to be eliminated. The corresponding UPF can thus be also A-dependent (like Z^(t + ), 
though sometimes this A-dependence occasionally disappears, like for Z~(i~( 2 ))). It is a partic- 

22 Indeed, if ni x ni matrix Li is identified in this way with the ni x ni matrix L2 (i.e. L2 has all the entries 
vanishing, except for those in the left upper m x m block), we have: tr„ 2Xn2 L2 = tr^xnjif for any L\ and k. 



ular integral representation (of the restriction of this UPF on the n-dimensional hypersurface) , 
arizing when n = N, that reproduces the original unitary matrix model per se. 

Usually in the unitary matrix lattice models of YM theories it is the weak-coupling (pertur- 
bative) phase is complicated. Our observation in this paper is that the strong-coupling limit 
in (at least some of) these models can be identified as a weak-field limit of certain versions 
of the GKM. Both these limits (strong coupled unitary matrices and weak field in the GKM) 
are reasonably simple. The outcome is that the interesting weak coupling limit in lattice YM- 
related models can be now identified with the strong-field limit of the GKM, which is rather 
well understood. 

We illustrated all these ideas with the simplest examples of the BGWM and IZ models, 
though the applications can appear to be much broader. Specifics of these models is that they 
are single-link theories and thus arise in the study of Yang- Mills theory only under restrictive 
assumptions (trivial topology in d = 2 or the mean- field approximation in higher dimensions). 
Another specifics - at least of the BGWM and even among the single-link models - is that it is 
essentially of the eigenvalue type. Going beyond the BGWM towards generic lattice Yang-Mills 
theories and c > 1 models requires broad generalization of Cartanian integr ability. 23 

The last thing to be mentioned is that the language of universal partition functions is the most 
adequate for consideration of "phase transitions" . Phase transition is nothing but a singularity 
of the UPF at some point (or finite-dimensional domain) in the infinite-dimensional space of 
parameters (= coupling-constants = time- variables). In a given physical system almost all the 
parameters (the form of Lagrangian) are fixed. Then if the line, associated with non-fixed 
parameter passes through the singularity, the free energy is actually singular at some point: the 
phase transition occurs. Singularity can be, however, easily avoided if some other parameters 
are allowed to change. 24 In this sense there are actually no different phases in the framework of 
UPF: there are rather just different asymptotics of a single analytical function. These become 
really separated phases only on the low-dimensional hyper surf aces, i.e. when most of external 
fields are switched off in a given physical system. 

These trivial remarks are important to keep in mind in application to particular models, 
which we discussed in the main text. The BGWM is known to have a phase transition in the 
N = oo limit. Actually this result was obtained in [3], when the space of parameters was severely 
restricted (J = a I) and the phase transition can look differently (or disappear) on other lines 
in the i-space. Moreover, unitary matrix integrals are, of course, usually finite for N < oo. 
This, however, does not imply that, say, Z^(t + ) has no singularities at finite values of iV: the 
statement is only that the n = N hypersurface tk = ^Tv(JJ^) h in the t-space avoids all these 
singularities (if any). Moreover, Z^(t + ) certainly has singularities as a function of N, which are 
again avoided by above hypersurfaces. For actual investigation of the phase transitions in the 
BGWM one can examine the large- N limit of Z^(t + ) by the method of ref.[10] (see section 4 of 
that paper or section 5.3 of [1]). Another thing to be taken into account is the switch between 
different UPF: from Z~^(t + ) in the strong-coupling "phase" to Z~(t~(~ 2 )) in the weak-coupling 
one. More detailed discussion of phase structures and singularities of r-functions is beyond the 
scope of the present paper. 

23 Cartanian r-functions (i.e. those of multicomponent Toda systems, including their KP, KdV etc reductions) 
are defined as correlators of the free 2d fermions and thus can be interpreted as determinants of 8 operators and 
described in terms of Grassmannian (Sato-Segal- Wilson theory). The necessary generalizations, first, substitute 
free fermions (fermionization of U(q) 1 Kac-Moody algebra) by generic 2d free-field theory (Wess-Zumino-Witten 
model), thus going from Segal- Wilson construction to generic description of determinant bundles over the universal 
module space. This should allow to go from "Cartanian" and thus essentially Abelian hierarchies to generic non- 
Abelian ones (i.e. from commuting Hamiltonian flows to those which from closed non-trivial algebra). Second, 
to go further beyond single-link theories (e.g. to those involving plaquettes), one should probably abandon 2 
dimensions and consider correlators in the free-field theories in d = 3. See [1] for more discussion of these issues. 

24 The most familiar example is of course elimination of the second-order phase transition by magnetic field. It 
is only at the line H = that the free energy is really singular. 
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Appendix 

This appendix contains proofs of several important formulas, which were omited from the main 
body of the paper, but can still be useful for illustrative purposes. 



In the first part of the Appendix we consider some subtle points connected with the transition 
from the external matrix variables A, J to the radial matrix variable M = JJ^ . This is essential 
for the derivation of the equation (3.5). 

Let us consider the following second derivatives 



04 &V 



kj 

which arise in the Ward identities satisfied by BGWM (see (3.4)). Now let these derivatives act 
to the function which depends only on M. Then, we can do the change of variables 
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for the derivative (A.l) and we used the definition of matrix derivative which is the derivative 
with respect to matrix elements of transponed matrix. Absolutely analogously one can do change 
of the variables in the derivative (A.l) 



d d _ d dM qs d d d d d d 
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Unlike the previous case this result can not be expressed in terms of radial matrix M. Nev- 
ertheless, it was just this expression which was used to derive Ward identity in the paper [6]. 
This is resolved by noting that the authors of [6] considered only this expression after changing 
variables to the eigenvalues A of matrix M. In this case, (A.3) can be really rewritten in terms 
of A and coincides with the formula (A. 2) expressed in terms of eigenvalues: 



This expresses the fact of completeness of these Ward identities which unambigously define 
unitary matrix integral (1.1) (this integral satisfies Ward identities with both possible derivatives 
(A.l) and (A.1)). 



In the second part of the Appendix we are going to discuss the leading term at the r.h.s. of 
(4.2). It can be evaluated as follows. Let 



- d 

K = tr$- — $XJ. (A.5) 
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The quantity to evaluate is: 
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Then 

1 1 n_1 ■ K m ■ 1 n_1 

k n = - : K n := S n , + -r V -tv^ n - m tr(^X) n - m = S n , + -r V fc m s n _ m . (A.7) 
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(A.8) 



Then 



or 



k(t) = 1 + /' k{t)s(t)dt, (A.9) 
Jo 



= k(t)s(t), 



k(t) 
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exp y s(t)tft = exp y — tr$ m tr($X) m = Det _1 (7 ®I- t$X ® $). 



(A.10) 
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